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I. INTRODUCTION 


A. BACKGROUND AND OBJECTIVE 

The purpose of this study™is to desceribe and anakyzewa 
numerscal summation of the Califommua coastalegone using ewe 
DIECAST (Dietrich-Center for Air Sea Technology) ocean model, 
earer™ 1 mproevedmaphyomes “Diet ri chyrrrooryy and seo" compare the 
observations/results with recent in-situ and remote sensing 
observational studies. The study is essentially a 
continuation of the study conducted by Akahoshi (2995) In 
Akahoshi’s study the DIECAST ocean model did not include a 
surface wind forcing but did include a simplified surface 
buoyancy forcing with damping to the  Levitus (1982) 
climatological summertime means of temperature and salinity 
and an equatorward-flowing surface jet at the northern 
boundary. Observationally realistic boundary conditions with 
realistic topography and coastal geometry were also 
incorporated into Akahoshis’ model study. In his paper 
Akahoshi concluded that the DIECAST ocean model simulation, 
utilizing realistic topography and coastal geometry with 
beendary «forciniqeealonerrm@canm™predice resus ==tnaewrarec 
consistent with observations to the first order. Akahoshi 
also conjectured that the inclusion of a surface wind forcing 


and ‘a8 weouthemme bouwmdary “forcing ine the form Gm a 


poleward-flowing undercurrent would likely bring model results 
more win limeewenwebservations. @ekahoshi also stated that all 
model forcings skewld inchudesanwemmual eyclesto peeduce the 
desired end result. Akahoshi’s findings were such that 
further study of ,the DIECAST ocean model seemed warranted, 
henee the ineorpomation of the impnewed medeimephysics. 

The DIECAST ocean model is a robust, state-of-the-art, 
high resolution, regional model. Primary asthengishns of athe 
model include its’ ability to generate and resolve fine eddy 
structure in three dimensions and eddy propagation, and its 
abilitym@eo be wmel@@ateds te any weean ocr@lake region with 
relative ease. The current simulation with the DIECAST ocean 
model incorporates significantly improved physics over that 
utilized by Akahoshi (1995). Significant improvements in the 
Simulation include surface wind forcing by the Hellerman and 
Rosenstein (1983) mean monthly climatological wind fields, an 
additional coastal wind enhancement (headland winds), surface 
buoyancy damping to the Levitus (1982) climatological monthly 
means of temperature and salinity, barotropic inflow and 
outflow at the northern and southern boundaries, and a 
baroclinic coastal jet at the northern boundary in phase with 
the surface winds. A more detailed description of the current 
version of the DIECAST ocean model can be found in chapter 


Lwien 


Be PREVIOUS STUDIES 

1 Enriquez and Friehe, 1995 

In their research, Enriquez and Friehe conducted an 
analytical and numerical study of the relative importance of 
wind@stressm(T) andthe cum! of@eherwinererness (Vx pon 
coastal upwelling/downwelling and its role in generating 
circulations near the coast, i.e. the time evolution of the 
mixed layer depth, A. The geographical location of their 
study was centered around Point Arena from 38° 12'N to 39° 12'N 
and from 124° 24'W to the coast. For their study Enriquez 
and Friehe used aircraft derived wind data from low-level (30 
m) flight tracks, and from vertical profiles north of Point 
Arena and south of the Russian River. The observed data, 
which showed enhanced wind stress and wind stress curl just 
offshore and south of Point Arena, was applied to a simple two 
layer, vertically integrated, linearized, reduced-gravity 
numerical model of coastal upwelling. The lower layer was 
dynamically inactive and the motion of the upper layer was 
represented by the first eee mode. Finite differencing 
was accomplished on a staggered Arakawa C grid using a 
leapfrog time advance scheme. No-slip conditions were applied 
to a 2 km x 2 km square grid coastline. Open boundary 
conditions existed on the northern, southern, and western 


boundaries. The analytical model was a simplified version of 


Ene wmm@merical model amd assumedya straight Coastline™ No-slip 
boundary conditions were applied at the coast and finite 
values of the variables were required far offshore. 

Rhinee. soluiuens @weae .obbaawncde frome bhey study =-— an 
analytical steady-sate solutionfewan analytical time-varying 
solution, and a. numericalesoliiticonmeRachwsolution showed that 
the application of wind stress curl near the coast had a 
Significant effect on the upwelling/downwelling regime. 
Specifically, curl not only increased the rate of upwelling 
and decreased the rate of downwelling but caused the 
horizontal extent of the effects of upwelling to expand beyond 
the euadeles application area. Additionally, loead 
(neighboring) extreme values of curl tended to merge into 
smoother large-scale.structures with less pronounced peaks. 
The magnitude of upwelling near the coast was uniform, 
decreasing farther offshore, and was strongest near areas of 
Mol i@lele wey eve oes: wie va maxima. 

Sustained high concentrations of positive wind stress 
curl observed near Point Arena are thought to be caused by the 
local coastal topography affecting the low-level wind flow. 
The study showed the importance of such concentrated wind 
stress to coastal upwelling. 

2. Burk and Thompson, 1996 


Burk and Thompson (1996) used this study to document the 


structure, and to investigate the mesoscale processes involved 
in the spring/summer time Low Level Jet (LLJ) that flows along 
the California coast from Cape Mendocino werPointCongepumen. 
A mesoscale atmospheric numerical model, Navy Operational 
Regional Atmospheric Prediction System (NORAPS), as described 
by Hodur (1987) with a modified physics package as described 
by* Burk and Thompson 4989), was Weed to condver@ehe™sure,. 
Optimum interpolation analysis was used for observational data 
assimilation. A control experiment was conducted by producing 
Sea Surface Temperature (SST) fields from the regional 
Expanded Ocean Thermal Structure (EOTS) and Optimum Thermal 
Interpolation System (OTIS) as described by Clancy et al. 
(1a9'9 0 jm 

In the@eestudy, Burk and Tihtempson found th@irehe iin Wad 
several regions of local maxima along the coast in the lee of 
capes and points. Associated with these regions of local LLJ 
maxima or “patches” are areas of enhanced upwelling and cold 
SST pools produced by the enhanced surface stress maxima 
caused by the topographically forced LLJ patches. 

Be Dorman et al., 1998 (In review) 

In their study, Dorman et al. CSOs ™ conetieted 
instrumented flightesurveys and real-time and retrospective 
numerical model forecasts of the California coast to determine 


the correct mesoscale structure of the Marine Boundary Layer 


(MET Ometesrne ocean, particularly in the vicinazty ofeeBeiunt 
Sur. Aircraft data was augmented with both fixed station and 
ocean sensing systems. iheeNavyes menhydrostatic Coupled 
Ocean/Atmosphere Mesoscale Prediction System (COAMPS) was used 
PbO provide, a Warger, scale, spatially .cOntinuous, three 
dimensional context of the marine boundary layer structure for 
the aircraft. 

One» relevant conclusion ~of «the Dorman etw alr. (1998) 
study, as it pertains to this study, Ys that the sebserved 
structure of the coastal MBL, supported by the model, is 
senonghys sugqgestavem=ewof wane imbounde=supererit hecaleeeMBL 
interacting with the topographic bend in the coast at Point 
Sur, forming a supercritical expansion fan in the lee (Dorman 
et wa lege 1 9SBix. Observationally, this correlates with their 
measurements of the lowest sea surface temperatures being off 
Point Sur near to, but not exactly under, the wind speed 
maximum. | 

4. Summary 

As described in the previous sections, the observational 
study conducted by Dorman et al. (1998) showed that the 
summertime influence of the subtropical high on the California 
coast was such that the equatorward flowing geostrophic wind 
developed a low level, concentrated, atmospheric jet on the 


leeward side of Point Sur. Dorman et al. (1998) theorized 


that such a feature may also Se present aby oulemet oes ie lee 
bends like Cape Mendocino and Point Arena. It is believed 
that the effect of this feature is such that the enhanced wind 
stress and its curl locally increase the off-shore Ekman 
transport and..pumping,...tlereby oincreasing gupwelling™ and 
possibly eddy generation and filament formation. As noted 
previously, the atmospheric model study conducted by Burk and 
Thompson (1996) supports the existence of these wind jet 
maxima at all capes and points along the California coast. 
Figure 1.0 is a representative forecast from their study of 
surface wind stress maxima in the lee of points and capes 
along the California coast. Figure 1.1 is the low level wind 
field around Point Arena from the study by Enriguez and Friehe 
(1995). The complete significance of such enhanced headland 
winds for the ocean is still unknown however, and further 
study needs to be conducted before definitive conclusions can 
be reached. The primary questions that need to be answered 
are; 1) Does the enhanced wind feature really exist at every 
promontory and cape? 2) What is the proper structure or shape 
of the headland wind jets? 3) What is the proper orientation 
of the jets? 4) What is their appropriate magnitude? 5) How 
do these features change seasonally, temporally, and/or 
diurnally? 6) How significant is this feature in the over all 


coastal dynamics? 


41N 





i : 3 0.1 3 ne 
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Figure 1.0 - Twelve-hour forecast of surface stress (N m‘%) 


var 1700 Bb 21 July 1997 (coneour Interval gs O.o2 N mm 
From Burk and Thompson (1996). 
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Figure 1.1 - Computed wind stress in pascals. Conteur 
interval = 0.05 Pa. From Enriquez and Friehe (1995). 


The present study, utilizing an idealized form of these 
strong coastal winds, is the first one to specifically address 


the question of their importance in forcing the coastal ocean. 
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Il. MODEL DESCRIPTION AND CONFIGURATION 


The following discussiom provaides. a description of the 
DIECAST ocean model and its configuration for the Simulation 
being studied. 

A. DIECAST MODEL 

The DIECAST ocean model is a primitive equation, z-level 
ocean/lake circulation model (Dietrich, 1997). eas 
hydrostatic, inceompressabley rigid=lmidiepaervaily Empigeit, 
end fully" conservative” The model is robust with real 
(unfiltered) topography, has very low (realistic) dissipation, 
and incorporates full thermodynamics. 

The model uses a pair of colocated grid structures for 
maximum computational accuracy and minimal numerical 
dispersion. An Arakawa A-grid is used for all computations of 
agwection, difiusion, fiction, @eGemiclis, ame bazcel mic 
pressure gradient forces. This greatly reduces numerical 
dispersion for the Coriolis term and a fourth order pressure 
Gradient cadcudeatiom amprovwes ives: accuracy omethe A—gagmd. 
Velocity is then interpolated by a fourth order method to an 
Arawaka C-grid structure where boundary velocity updates, 
incompressibility constraint calculations and removal of the 
Vererea lily integrated divergence are accurately and 


efficiently performed. Velocity values are then interpolated 


dete 


back to the Arawaka A-grid. This improvement to the DIECAST 
©Oceenpmete mma conpamed LO the study by Akahesia (1995), results 
in far less numerical dispersion with a much improved handling 
of eddies and fine scale gradients. 
iB: MODEL CONFIGURATION 

a Domain 

The domain of the model is from 32.0°N to 42.0°N and from 
132.5°W to the coast. To remove possible boundary influences 
on the analyzed fields, the analysis domain has been confined 
mo themrequons from3495eN towm41. 0°NS@and @iromel2650°W temEne 
SCelieeornia cogs. Sts Sxtends Mrom@Pennt Conceperen to Gust 
noren of Cape Mendocino end irene Ene Caliiounta coast Co aver 
600 km offshore. The coastline topography includes the major 
headlands of Cape Mendocino, Point Arena, Point Reyes and 
Point Sur (see Figure 2.0). These headlands will become an 
mimoctant Geplc of discussion below. 

2. Resolution 

The longitudinal/horizontal resolution of the model is 
locally fixed at = longitude. Latitudinal resolution is 
variable such that the grid spacings are equal, Ax = Ay, where 


Ay = phi eco s (Gia: Ths approximately equates’ to 


eZ 
hx SPA C2 SP at. 34 .S°NWand AY SyreanG. 98a? at 4ae" The 
model has 20 vertical layers of resolution with the uppermost 


layer being 21 m thick and expanding down to a 689 m bottom 
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Figure 2.0 - Analysis domain depicting model topography levels 


10 through 20. Levels one through nine were omitted to reduce 
Slultter. ~Wepbaris in mMerens. 


irs) 


layer thickness (see Table 2.0). 
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Table 2.0 - Model layers with corresponding depth in meters. 
Az is the layer thickness in meters. 


Heghemmpe-coLULHONnS ins tne Uppenwocean,. 4.¢. ssme Beem 
Vomuicclimecemenwditterences, were designed ~to resolve. —ac 
lanje uy emercal Grdglchits present in eEhe upper ocean: In 
general, the high vertical resolution provides for a more 
realistic representation of topography (see Figure 2.0). The 
coastline and bottom depths were derived from the ETOPO-5 
topographic data sets (NOAA, 1986). 

Se Surface Forcing 

a. Climatological Surface Winds 
The climatological surface wind stress forcing comes 


from the Hellerman and Rosenstein (1983) data base, which 
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gives monthly averaged stress components. These are 
interpolated to the modelies. guid. using Same ei acia ene 
interpolation scheme. The monthly values are interpolated in 
time to daily values and then modified by an idealized 
headland wind enhancement (below) to define the total wind 
hOrCiINgG at thessumprace, 

b. Headland Winds 

In the present simulation, an idealized enhancement 
of the surface wind stress in the vicinity o£ the coastal 
headlands is added to the climatological wind stress to 
provide the total wind stress forcing in the simulation. 

To do this, a headland wind enhancement factor is 
constructed by accumulating an idealized representation of the 
individual effect of each of four headlands (Cape Mendocino, 
Point Arena, Point Reyes, and Point Sur) om the surface 
stress. The resulting model wind stress is given by 

U aie (1 AE) (x,y ee, 


Vneve=L. Lsethe monthly mean climabologucal vameh Sieicessmet 


Hellerman and Rosenstein (1983) as described in the previous 
section, A(t) is the time-varying normalized amplitude of the 
additional headland wind effect and F(x,y) 1s its spatial 
pattern. Thus, the headland wind effect is a simple time- and 
space-dependent scalar enhancement of the climatological wind 


stress. This enhancement therefore does not change the wind 


i 


stress direction, which remains the same Sls. 1 O Mavehawnak 


1s currently known about the observed headland wind phenomena 
fogeme Cesc, ., 1220), Ee normalized amplitude function A(t) 
is specified to be non-zero only during the warm season of the 
year, March through September with a peak in May/June, the 
time of the strongest equatorward winds (Figure 2.1la). Also 
following observations, the spatial amplitude function F(x, y) 
is modeled using a Gaussian-like function centered just south 
of each headland. The spatial spreading and orientation of 
the fan-like wind patch 1s controlled by varying the length 
scales downstream of the headland. The enhancement function 
F(x,y), which describes the cumulative effect of all four 
headlands, is shown in figure 2.1b where the individual 
contribution from each of the four headlands is clearly seen. 
Since Point Arena and Point Reyes are sufficiently close to 
one another, the effects of these two headlands combine to 
give the greatest total wind enhancement off San Francisco 


Bay. The model wind stress in May/June is therefore up to 2.3 


pane. n thatWeream Figure 2.2 shews whe resultwng total 


surface wind stress fields at four representative times of the 


meat. 
Ce Surface Buoyancy Damping 


Surfacesbuowency seReimge!s veuy much simpdyetied and 
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Annual Cycle of Headland Wind Stress Enhancement Factor 
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Figure 2.1 - a) Monthly change in amplitude of headland wind 
stress factor. b) Magnitude and structure of headland wind 
stress factor. Non-dimensional units. 
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Gemsists of dampangy the surface fields. of temperature and 
salinity back to the annual cycle of climatological Levitus 
(1982) values on a time scale of 60 days. McCreary et al. 
(1991) observed that doubling the damping time scale from 40 
to 80 days in his model study produced similar results, while 
only taking longer for his model to reach equilibrium. Thus, 
the exact value of the damping time scale used is not 
Gencidercds tosbescritical, toe the sods 16m: 

Numerical damping of all fields is accomplished with a 
simplesdiftfiusion coefficient of 10° cm*/s which ineseascsmiag 
90% in a weak sponge zone at the western boundary. Tangential 
velocities near the boundary are also damped. Additionally, 
a weak biharmonic filter is applied to the barotropic velocity 
and the density field in order to suppress 2Ax waves. 

4. Lateral Forcing and Boundary Conditions 

All boundaries are partially open with the computed 
normal boundary velocity (NBV) determining inflow and outflow. 
On the western boundary, outflow of Rossby waves is enhanced 
by adding an additional iL cm/sec apparent outflow when 
determining the inflow/outflow condition. In all cases 
temperature and salinity are advected in or out, and momentum 
is damped as noted above. Initially, NBV is the specified 
normal boundary velocity, but during the simulation NBV is 


continually nudged towards the annual cycle of specified 
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February Wind Stress Contours and Magnitude Vectors 
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c) 


August Wind Stress Contours and Magnitude Vectors 
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Figure 2.2 - Surface wind stress contours and vectors. 
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May Wind Stress Contours and Magnitude Vectors 
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d) 


November Wind Stress Contours and Magnitude Vectors 
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mel@eities Onma 15 day time scale. The specified boundary 
Values Ot mommal velocity consist of several parts. The basic 
VelOcicy Jseee- NOrmalsgcostro phic baroclinic cuxnsent comemited 
from the anmmgaal cycle of temperature and salinity on the 
boundaryegiven by Levitus (#982). 

A specified barotrepic flow ws addé@ to the above 
geestresnic velocities at oth the northern aad southern 
boundaries. The barotropic normal velocities are composed of 
smoothly fitted values taken from the general circulation 
model of Semtner and Chervin (1992). The Semtner and Chervin 
model shows an annual cycle of barotropic flow through the 
California Current, northward in the summer and southward in 
the winter (Figure 223). ~In this study, ie across—-boundary 
barotropic transport is confined to the continental slope and 
modeled as a sinusoidal oscillation over the annual cycle with 
an amplitude of 5 Sv and a mean of zero. 

In the surface layer, the wind Be aven Ekman flow is 
allowed full freedom to change outflow. For subsurface layers 
the western and northern boundaries are constrained so that 
the mean value for each layer varies with the annual cycle 
boundary means. Each subsurface layer is also constrained to 
have no net divergence. 

Dp. Coastal Jet Forcing 


At the northern boundary, the Levitus (1982) 
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Annual Cycle of Poleward Barotropic Transport 
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Figure 2.3 - Annual cycle of barotropic flow through the 
Saldfornia Current. Normalizec™@cmrve"ol appl reamveluest 


Za 


See eo ea ecm pemaeure andwysalinity forcing, and the 
EGSU4bing Gee@strophic velocities described above, are modified 
to include an observationally realistic equatorward-flowing 
upper ocean jet. Ths» 1s edeeenplished bywadjusting the 
ViemEtcal Sseceions Of Mevitus temperature and salinity at the 
bountden,. brquie 284 “shows the “resm@lting “adjusted 
climatological temperature, salinity and v component of 
velocity at 42.0°N averaged over the month of May, the period 
of peak intensity and the month of November, an opposing month 
of relatively weak intensity. The jet (Figure 2.4c) is 
structured as a Gaussian jet with a core velocity of 30 cm/s, 
a horizontal scale of about 95 km, and a vertical scale Be 
abeutes 750° m. It 1S positioned adjacent to the continental 
slope in about 3000 m of water. The jet structure is intended 
to approximate the coastal jet observed in June 1987 during 
the Coastal Transition Zone (CTZ) experiment (Kosro et al., 
1991). The time variation of the coastal jet (Figure 2.5) is 
prescribed to metch the annual cyche of the baroclinic #et in 
the California Current as found in the Semtner and Chervin 
(1992) simulation that was also forced by the Hellerman and 
Rosenstein (1983) wind stress. The coastal jet is prescribed 
simply for consistency at the boundary, since we expect such 
a coastal jet to be generated immediately inside the domain in 


response to the wind forcing. 
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May Vertical Temperature Profile at 42 Degrees North May Vertical Salinity Profite at 42 Degrees North 
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Figure 2.4(a-c) - Boundary conditions at 42° north for the 


month of May. a temperatume s@@)r, bh) salimteye(psu)y chev 
welocity componemt (em/s) . 
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November Vertical Temperature Profile at 42 Degrees North 
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Figure 2.4(d-f£) - Boundary conditions at 42° north 
month of November. a) temperature (°C), b) salinity 
c) v velocity component (cm/s). 
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November Vertical Salinity Profile at 42 Degrees North 
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Figure 2.5 - Time variation of baroclinic jet at the northern 
boundary. Normalized curve of applied values. 


Zo 


oe Time frame 

The model was run for a total time period of six model 
years. Initial conditions consisted of one degree resolution 
climatological monthly means of temperature and salinity from 
Levitus (1982). This data was interpolated and gridded to our 
models matrix to provide a full field of T,S values for the 
initial January O01 start time. As described above, the 
monthly values of T,S on the lateral boundaries (adjusted to 
produce the coastal jet at the northern boundary) were then 
used to force the annual cycle. In addition, the surface 
values were used to provide the annual cycle for surface 
BeSteraelon. 

The first year of the model run was considered to be the 
spin-up period while years two through six were considered to 
be in neéar—-equiliybrium. Figure 2.6, a time series plot of tae 
daily mean eddy kinetic energy at the surface, clearly shows 
the rapid adjustment period during year one and the slightly 
downward sloping trend in energy levels from years two through 
six. Because of this trend the model was considered in near- 
equilibrium durdipe years two throdenh siganed@this consiiveume 


the majority of the analyzed data fields. 
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Figure 2.6 - Daily mean Eddy Kinetic Energy at the surface. 
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III. SURFACE VARIABILITY 


This chapter will analyze the model simulated surface 
variability and compare the results to observations. 
specifically, the characteristics of dynamic sea surface 
height and eddy kinetic energy will be studied and compared to 
a recent comprehensive observational study conducted by Kelly 
et al. (1998). 

A. DYNAMIC SEA SURFACE HEIGHT 

The resulting five year mean dynamic sea surface height 
(SSH) fromthe “Simulation (Figure 3.0) shows equatcorward 
geostrophic flow throughout the region, with the broad 
California Current well offshore and the largest SSH gradients 
near the coast from north of Point Arena to south of San 
PranewmescO Bay. the “Sea field has a Maine Sielies oF 
approximately -O.izZ m just north of Cape Mendocino near 41°N, 
124.5°W and a maximum value of 0.12 m near 34.6°N, 128°W giving 
a total SSH difference of 0.24 m. Kelly et al. (1998) 
observed a SSH minimum near 36°N, 123°W (see Figure 3.0) with 
a SSH difference of 0.22 m over the same domain. 

A significant difference in the model study compared to 
observations is the tight gradient of SSH contours that run 
parallel, “and close to the eeasc. The -O-12 m corcour Eins 


adjacent to the coast from the northern boundary down to 
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Mean Dynamic Sea Surface Height (m) 


40 


39 





=e 
= 
oO 
z 
g 38 
mo 
® 
Ss 
® 
TC 
= 
® 
nay 

36 

35 

-128 -127 -126 -125 -124 -123 -122 -121 
Longitude (Degrees West) 

Pigure’ 3.0 - Mean Dynamic Sea Surface Height, in meters: 


Genweur intermeal is 0.02 mee The “X%enea@r 36°N, 126°W marie Ehe 
location of minimum sea surface height observed by Kelly et 


al. , giclee. 


30 


37.3°N, just south of the San Erancisco Bay areamandetne eam 
m contour runs parallel tomwit from about 9399S -N to sae 
36.3°N, just west “Gm Ppowmememr= Tre -0.08 m conisgur becuase 
its run near 39.4°N and goes southeast down the coast to the 
eastern boundary. There is a large trough, or sea surface 
depressiom, extending from Cape Mendocino southwest toWabout 
39.5°N at the western boundary. A smaller rise and depression 
pair appear southwest of Monterey Bay. This pattern is in 
eontrast to Kell s findwngs for the mean SSH in whidema 
Single minimum value is located along the coast just west of 
Monterey Bay with rather uniformly spaced SSH contours 
generally “parallel to the coast but archimgeecyclom calmly 
around the minimum. The deep trough in the simulated SSH 
field south of the Mendocino escarpment is simply not seen in 
the observed field. 

Comparison of the model mean dynamic sea surface height 
(Figure 3.0).wathe thesemedel mean sea sumiace pempematume (SST) 
field (Figure 3.1) shows a high positive visual correlation 
between SSH and temperature, especially near the coast and 
west of the Monterey Bay and Point Sur regions. Particular 
patterns to observe are the troughing south of Cape Mendocino 
and west of Point Sur and the ridging west of Point Arena and 
Monterey Bay. The minimum Sodeeais located along Che™coast, 


coinciding with the minimum SSH as expected with coastal 
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Figure 3.1 - Mean Sea Surface Temperature in 
iiwerval is a°C. 
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upwelling. Seasonal comparisons of sea surface height and 
temperature will be conducted in the next section. 

The strongest mean geostrophic currents, computed from 
the mean SSH field, run adjacent to the coast from about 
39.25°N to about 37.0°N reaching a maximum velocity of about 
0.40 m/s near 37925.) 122. 75°Weethes Compares to Kelly et 3am 
(1998) who observed a predominately southeastward one year 
mean geostrophic jet, computed from drifter and hydrographic 
data, centered along approximately 36°N, 125°W to roughly 31°N, 
122°W, with a maximum southeastward velocity of about 0.13 m/s 
at@sem, 1265. At 37°N; i25°W, the meeaimum observed mean 
geostrophic velocity is southwestward at about 0.09 m/s. 
Thus, the model simulated mean surface current maximum of 0.40 
m/s is at least three times larger than the maximum observed 
mean value computed from combined surface and drifter data and 
elimietological hydrographic data: 

a, Combined Seasonal Dynamic Height and Temperature 

For this study the seasons are defined as follows: Winter 
(January, February, March); Spring (April, May, June); Summer 
(July, August, September); Fall (October, November, December) . 
The following discussion, referring to Figure 3.2a-d, shows 
the colored temperature field overlain with contoured dynamic 
sea surface height for each season. 


Starting with the winter season (Figure 3.2a) one sees 
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Winter Temperature and Pressure Fields 
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Figure 3.2a - Winter temperature field (colored) in °C with 


the winter surface dynamic height field in cm overlaid. 
Contour interval 1S 2 cm. 
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temperatures in the range of about 11°-12°C at the coast, from 
the northern boundary west of Cape Mendocino to south of Point 
Sur, 6 about I73@iimear 35 Nee? Sow. Streng rigging in the 
temperature ander essure Sat awees evident orle@m@matang at 
about 36.5°N and extending toward the east and northeast. 
Local dynamic height highs are located about 140 km west of 
Point Arena and Monterey Bay. @roughing in the temperature 
and pressure fields is also evident extending southwest from 
Cape Mendocino and from Monterey Bay/Point Sur. A small 
trough is present immediately west of Monterey Bay. A closed 
low/cyclone is present near 39.5°N, 126.5°W. Right along the 
coast the dynamic height and temperature fields exhibit a 
relatively relaxed, weak pattern showing strength only near 
Point Arena and just south of the San Francisco Bay region. 
Equatorward flow is evident far offshore in the southwest 
quadrant and adjacent to the coast. By comparison, Strub and 
James (1998, in review), who analyzed SSH fields from four 
years of altimeter data, observed a dynamic low in the height 
pattern west of Point Arena, without the ridge of high 
pressure seen there in this study, and a low near 35°N, 123°W, 
Seighit ly “tow rhne=seuth Gf the Point ‘Sur trough in’ thitsmstudy . 
Strub and James observed a total temperature range from about 
10°C to 15°C, smaller than observed in this study, but with a 


Simi lastedilst rat ilem Their height gradient was also less 


a0 


Cid S Stuay, i2 sem verses 22 cm, wath osimilar 
distribution as discussed as above. 

Durie thes spring seasom (Frgure 3:25) both the 
temperature and the dynamic SSH gradients reach their maximum 
values - a 7°C southwest SST gradient with a 32 cm southwest 
dynamic SSH gradient. The coldest upwelled coastal waters 
only go as far south as about 37°N, increasing in temperature 
south of that. The cold plume west of Cape Mendocino in the 
winter has migrated to the coast in spring. The gradient of 
the western boundary SSH ridge has relaxed significantly 
however the closed high/anticyclone near 37°N, 124W, present 
in the winter season, remains and has intensified by at least 
6 cm. The gradient of the troughs, also present in the winter 
season, extending from Cape Mendocino and from Monterey Bay 
have increased and a closed high has developed near 35.5°N, 
122.5°W from a weak ridge present in the winter season. 
Upwelling along the coast from Cape Mendocino to Monterey Bay 
has» reached its iaeRimum prmeeticing@e the “rightest, “Wee 
concentrated SSHegradient of the four s@a@sons. Seasonal mean 
geostrophic velocities along the coast have also reached their 
IMaximums of about 36 cm/s just west of Cape Mendocino to about 
50 cm/s between Point Arena and the Monterey Bay region. 
Equatorward flow is more uniform, paralleling the coast near 


shore from the northern boundary to about 37°N and far off 
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Figure 3.2b - Spring temperature field (colored) in °C with 
the spring surface dynamic height field in cm overlaid. 
Contour interval is 2 cm. 
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shore south of 37°N. Again, Strub and James (1998, in review) 
observed a mearly ideng@iecaleemetribution pattern of SST and 
dynamic SSH, ditiftering from the simulations in the magnitude 
of the observed gradients - southwest gradients of temperature 
(about 6°C) and dynamic height (about 20 om). An 
approximately 2 Nem =eceestroplmemsl|ow parallels the entire 
coast .@eegain, ~thewmedel gradivenesmee SSH and the resulting 
mean geostrophic currents are larger than observed values. 
The SUMMGT pattemiaowet soot ciemewes (igure 3.2c) have 
become significantly relaxed compared to the spring patterns. 
From the coast to well offshore the SST ranges from about 
9°-15°C, and the SSH ranges=rron —lwemetom=emem-:  Upwelling 
along the coast has become less pronounced though its extent 
along the coast remains relatively constant. aige ridge— 
trough-ridge pattern west oof Point Sur has relaxed 
considerably from the spring season, suggesting that the 
strongly seasonal headland wind jet may be the major driving 
factor. Strub and James observed a similar relaxation of the 
temperature and height fields but with the coastal jet 
migrating farther offshore and becoming more diffuse, rather 
than just weakening in place as seen in the simulation. The 
observed extent of coastal upwelling has significantly 
decreased by summer and the coldest temperatures in the 


upwelling region have slightly increased. The observed 
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Figure 3.2c - Summer temperature field (colored) in °C with 


the summer surface dynamic height field in cm overlaid. 
Contour interval is 2 cm. 
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temperature range across the domain were about 11°-19°%, 
somewhat warmer than that found in this model simulation. The 
mean SSH shows an observed range of about 18 cm in summer, 
much lower than the 24 cm SSH Bagge found in this simulation. 
The observed fields in the fall season (Figure 3.2d) have 
become completely relaxed showing much less structure and 
detail than in the previous seasons. Upwelling along the 
Geast has reached 7is minimum extent ;Wemusting only in a small 
band from Point Arena to Point Sur. The isotherms are 
distinctly oriented northwest/southeast and range from about 
11°C near the coast to 17°C in the southwest. Dynamic SSH 
ranges from -10 cm along the coast to 12 cm near 34.6°N, 
128 Came The ridge west of the Monterey Bay remains 
persistent. The trough south of the Mendocino escarpment near 
39.5°N, 126.0°W has closed off and is interacting with a weak 
coastal ridge slightly off shore adjacent to Point Arena. 
Strub and James show a Similar pattern of relaxation with 
nearly parallel alignment of the temperature and height fields 
with the coast. However, the observations show that the 
coastal jet has weakened and moved farther off shore, 
resulting in a broader coastal area of cooler temperatures 
with northward flow developing along the coast from Point 
Brena to 4Ul0°-Ne ine model belavilor, in whitch tne so Ulm crae 


coastal current simply weakens in the fall, is thus different 
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Figure 3.2d - Fall temperature field (colored) in °C with the 
fall. surface dynamic height field in cm overlaid. Contour 
interval 1s 2 cm. 
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from the @@bservedmbehavior in which We Goeasts ema me 
propagates offshore. Observed temperature and dynamic SSH 
ranges in Strub and James are about 7°C and @20en 
respectively, comparable to that of this simulation. 

eee Discussion of findings 

One difference between our findings and those of Strub 
and James (1998, in review) is that the model fails to show 
the seasonatepaeeemimor Ememcoastal jet moving off Sshiore ame 
relaxing from its strongest season, spring, to its weakest 
season, winter. The model shows a strong coastal jet, 
essentially from Point Arena to MontereyeBay swireheenersommung 
and stmmem se@@sonsm being its sttong@est,. relaxing ame 
Gepeeacting through the fall, omni to begin Sstrengesenmnc game 
Spanding again in the beginning) of sprang. ~DBhileaimayeseus 
function of the strength and duration of coastal upwelling 
caused by the idealized headland wind forcing. As noted 
earlier, this is an area of isin uncertainty at the 
present time. 

The pattern of substantial ridges and troughs along the 
coast in the simulation is not so evident in the observational 
study of Strub and James (1998) where meanders in the along 
shore flow are much weaker. This may also be a characteristic 
of the model producing overly strong sea surface temperature 


and dynamic height gradients causing the resulting strength in 
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the ge@strophic flows to interact lwithmehesbecromutepocm ih, 
thereby enhancingmany natural cemdeneyemeemmidge Ore nouan 
(Refer to Pigtire 2% 0)™™ One othe poSsmevlity Pemmehe Cmigemeed 
features may be artifimenal inttwences from deem ortshere 
boundary. These. possible boundarye -influiemeées wilwee be 
discussed further in the next section. 
B. EDDY KINETIC ENERGY 

so Mean Kinetic Energy 

The “tearrecayarelimecic energy fieTd (ENB etrigure 3. 
shows two main regions of relatively high energy at the 
surface: Region one, located south of Point Arena near 38°N, 
123.5°W has a local maximum value of 0.03 (m/s)* with a sub- 
megtOm= valle Gt OPO25 im si about degree latitude to its 
Seutch; @Region two,@located south™of "Peint Sur near 35.5, 
122.5°W has a local maximum value of 0.035 (m/s)*. Minimum 
values of 0.005 (m/s)* are located at the northern boundary 
and along the immediate vicinity of the entire coast. This is 
in contrast to Kelly et al. (1998) who observed a single local 
maximum value of about 0.045 (m/s)* farther offshore near 
TS°N , 225 SW ater keayen Figuree3 9S) with a ningun" emtour 
of 0.015 (m/s)* along the western boundary from 127.5°W to 
128°W. The EKE values shown in Kelly et al. (1998) are based 
on three years of surface drifter data which unfortunately 


only sampled the region seaward of 124°W. The EKE 
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Figure 3.3 - Time mean Eddy Kinetic Energy in (m/s)*. Contour 
mitervel Wee0,005e(m4s)*. “The “x” “at approx. 36.5°N, 12695°%H 
marks the location of the maximum value observed by Kelly et 
ay Aes». 
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observations therefore do net include the Goastal areaues- mer 
124°W where the model simulated EKE is largest. 

The simulated EKE gfield in @Pigure 3 .3 johows jamie weno 
minimum near 125°W at most latitudes. West of about 125.5, 
and south of about 40°N the EKE increases westward, showing 
thatemenergy %1s ,ea1ther being Ge@meratcdm locally om is 
trapped/not efficiently dissipated at the western boundary. 
The normal expectation is that energy generated at or near the 
coast would propagate westward and either dissipate and/or be 
dispersed through a deeper column of water leaving lower 
enemey levels at the western boundamy gion amy Gaven digger seit 
is also expected that EKE would only be generated at or near 
the coast as a result of coastal upwelling and wind driven 
eddy generating dynamics at the surface. me. EtEurther 
investigate the seemingly abnormally high energy levels near 
the western boundary the EKE at two lower levels were 
examined. Figure 3.4 not only shows that higher energy levels 
exist west of about 124.5°W at the deep levels of 9 (392 m) 
and 14 (113 m), but energy levels are actually increasing 
toward the western boundary. This is not consistent with the 
normal expectation that energy levels should be decreasing 
westward. Thus, a working hypothesis at this stage is that 
the westward increase in EKE west of about 125°W at all levels 


is due to an artificial reflection of energy from an 
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Mean Level 9 EKE (m/s)? Mean Level 14 EKE (m/s)? 
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Figure 3.4 - a) Mean Eddy Kinetic Energy at level 9 (362 m). Contour interval is 0.0002 
(m/s)*. Bb) Meam EKE at levelel4 3 ial aie Contome: Bntervak is O@9001 Gi/s)-. The 
innermost dashed lines adjacent to the coast represent the topographic demarcation at that 
level. 
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incompletely “open” western boundary. 

aa Temporal Changes In Eddy Kinetic Energy 

To avoid possible influences from the western boundary, 
the EKE fields were further analyzed in a coastline relative 
frame of reference. This “domatinm@of analysis extends se) 5 
degrees of longitude (about 330 km) westward from the 
coastline at each latitude. This region extends well into the 
deep water beyond the continental slope (Figure 2.0), but it 
excludes the region of artificial western boundary influence. 
Because the continental slope extends westward with increasing 
depth, the easternmost data point for lower levels (the 
coastline of those levels) also moves westward. This means 
that the horizontal width of the domain at each successively 
lower level decreases slightly as the western edge of the 
domain is fixed (330 km west of the coastline at the surface). 

Figure 3.5a shows the time evolution of the EKE at the 
surface (level 1 (10 m)), averaged over this coastline 
relative domain. The time series shows a consistent pattern 
of peak seasonal energy levels occurring in the month of 
May/June. The maximum observed value is 0.0241 (y's) -. 
Maximum values occur in late spring/early summer and appear to 
coincide with the annual cycle of headland winds (Figure 
2. ey The minimum value of 0.0069 (m/s) * Sceurs = ii late 


fall/early winter and also appear to correspond to the annual 
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Daily Mean EKE (Level 1) 
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cycle of the headland winds. As noted in a previous section, 
the EKE time variability appears to be seasonally adjusted at 
about 0.0133 (m/s)* for an annual mean. These values 
correspond to eddy current speeds of about 8-16 cm/s with a 
mean of about 12 (cm/s). This compares with the two years of 
BRE from drifter “data in Kelly et ail (1906) wites sitars 
maximum valves Sof abowe) 0.02 (m¥S)* Seccurmimey in late 
summer/fall and minimum values of about 0.01 (m/s)* occurring 
in the early spring. As noted earlier however, the EKE values 
in Kelly et al. (1998) are averaged over a 700 x 900 km region 
entirely west of 124°W (i.e. well seaward of the slope) 
whereas the EKE values in Figure 3.5 are averaged over a 
region within 330 km of the coast. The simulated EKE reaches 
its maximum value earlier in the year (spring) than the 
observed EKE (summer/fall) because of the east-west difference 
in averaging domains (Simulated verses observed) and because 
the EKE has a pronounced westward propagation as described in 
the next section. 

Figure 3.5b shows the time evolution of the EKE below the 
main thermocline (level 9 (392 m)). As expected, the seasonal 
Signal at the surface appears to be reflected in the lower 
level at reduced amplitude and time lagged. The peak energy 


level now appears to be in the late summer/early to mid fall. 
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Daily Mean EKE (Level 9) 
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Sic Propagation of Eddy Kinetic Energy 

One of the most robust features of Ches@EKE anaiiysi caer 
Kelly et al. (1998) is the distimct westward propacaci on me 
decrease of EKE away from itS source near the coast. ne 
examine this characteristic 1nMehemmedome canteen Fi Gime 
3.6a shows a time distance (longitude) plot of the along-coast 
averaged EKE abltuhe Surtaces “Sonor odeaco le mmene analysis is 
performed over the coast relative domain. The left edge of 
the plots in Figure 3.6 are 330 km from the coastline (at the 
surface) while the right edge of the plots are at the 
coastline appropriate for the given level. The first evident 
feature is that high energy bands occur on a periodic basis, 
about every June, as expected. This EKE generation is clearly 
in connection with the spring/summer upwelling regime produced 
by the surface wind forcing. A tendency for offshore 
propagation with time is also evident. The peak EKE band 
oGeurs between about 35-185 km west o@ shore, with) max liam 
values occurring between 50-150 km west of the coastline, 
consistent with the annual mean EKE field of Figure 3.3. The 
westward propagation rate is estimated to be on the order of 
myo (2). KM er Cay OL eooUe 0.027 m/s. oye coMmpear 1som,, Kerry 
Se alin (192e" mi seTeee that peak EKE values at the surface 
also occurred in the June/July time frame near the coast and 


the westward propagation rate was estimated at 0.03 m/s, 
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Figure 3.6a - Time and distance evolution of Eddy Kinetic 


Energy at the surface (level 1 (10 I) i) eee bens 9s Pom bec eam 
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somewhat faster than observed in this study. These speeds of 
propagation are consistent with the speed of non-dispersive 
first mode Rossby waves at this latitude. 

Time distance plots Of ERE @ae levels so (3925m) ano 
(1139. m), SHGwn Gn bigqguresms. Gillie ae eo continue wee. Show tne 
westward propagation of EKE as expected. Most importantly the 
largest energy signals occur progressively farther offshore 
with depth, consistent with the energy moving out from the 
coast and penetrating downward with time. This vertical 
redistribution of EKE to the deep ocean west of its source 
near the coast partly explains the westward decrease of EKE at 
the surface that was documented in the California Current 
region by Kelly @@ al. (1998) and Strub amd Uaneemos). Mie 
strong seasonal signal apparent at the surface becomes 
incoherent and gets blurred with depth, where the time 
interval between energetic events appears to be of the order 


Ge enly i-2 mompns. 
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Figure 3.6b - Time and distance evolution of Eddy Kinetic 


Energy at level 9 (392 m)). EKE is plotted in days verses km 
from surface coastline. EKE is measured in (m/s)*. 
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Energy at level 14 (1139 m). EKE is plotted in days verses km 
from surface coastline. EKE is measured in (m/s)?*. 
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IV. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 


An updated version of the DIECAST ocean model was used to 
conduct a numerical simulation study of mesoscale variability 
in the California coastal zone. Improvements to the model 
included an annual cycle of surface wind forcing, headland 
wind enhancements, surface buoyancy damping, and time phased 
baretropac and” bameclimic acYessmbourdary meme - Speeaitic 
model results were then compared to recent in-situ and remote 
sensing studies. 

A. SUMMARY OF RESULTS 

The observed mean dynamic sea surface height generally 
shows southwest geostrophic flow north of Point Arena and 
southeast flow south of Point Arena. Model simulated height 
values run generally parallel to the coast, but with a trough 
immediately south of Cape Mendocino, contrary to observations. 
The total height range over the entire model domain is 0.24 nm, 
about 10% higher than observations. The strongest annual mean 
geostrophic currents are about three times stronger in the 
Simulation than observed. There is also considerably more 
structure and along shorewvariability im the simulated mean 
dynamve™" nervenet a1 dq thanmineehewebserved height” field. 

The seasonal sea surface height and sea _ surface 


temperature fields produced by the model show considerable 


of 


Sill iar ios tn oermceure tO Observation Nowever, the model 
Shows stronger gradients in both the temperature and height 
melds throwe@howes the four s@asons than®is observed. The 
southward coastal jet, always present in the model simulation, 
remains quasi-stationary and close to the coast. The jet 
pulses in Stréngthy dwem Stremgthens:' andm weakens @ever the 
Seasonse vise relaxiag and Smevangmoff shores during §the non 
upwelling season and then reforming adjacent to the coast the 
next spring. 

The simulated mean eddy kinetic energy field shows 
considerable structure and variabil] iseyjpeewart het wosmetin regivens 
Ommeccdy Veiemgy PwecuiclLionmmadjacent to the coast, one jUSe 
south of Point Arena and one south of Point Sur. Near shore 
observations of EKE are not available for comparison but an 
offshore study (Kelly et al., 1998) showed a single maximum of 
eddy energy west of Monterey Bay, in between the two regions 
the model simulated eddy kinetic energy maxima. Eddy energy 
decreases between about 124°W and 126°W then re-intensifies 
west of 126°W to the western boundary. The investigation of 
this phenomena at lower levels revealed that there is likely 
to be western boundary influences contaminating the simulated 
data from the western boundary eastward to about 126°W. 

Because of this, a coastline relative domain was 


developed to minimize the possible influence of the western 
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boundary on the analysis and interpretations. 

The temporal changes in eddy kinetic energy averaged over 
theseoastal region shoWéd”"a demmearee strong seasonal = ceyel——im 
phase with the coastal winds and in agreement with the 
analysiswot Kellyseteaiv (1998) Thessamulated energy Vevels 
at the surface are somewhat higher than observations but the 
differences are not considered significant since the domains 
bemmg comparedearessemewhated»iferemt. efhe saimulbeved phase in 
amplitude appears to coincide with the surface wind forcing. 
The surface energy signal appears to be reflected at lower 
levels, weakening and lagging in time as expected. 

The tameslongytude “evolutions of ~eddy ssinebiem etergy 
averaged along the coast shows a seasonal signal and westward 
propagation consistent with bothwtemporalychanges inweddy 
kinetic energy and westward propagation seen in observations. 
Our westward propagation rate is estimated to be somewhat 
slower than in observations, but the difference is not 
considered significant. The downward spreading of eddy energy 
is evident in the simulation with peak EKE levels occurring 
further and further offshore with depth, while the seasonal 
Signal of EKE becomes indistinguishable at depth. 

B. CONCLUSIONS AND RECOMMENDATIONS 
Results indicate that the use of improved model physics 


can produce ocean variability generally consistent with 
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Seeer 2 erons Vineet Cali fornta coastal zone, but “With 
important differences as well. 

The nesult tat the model»preducedwa stronger gradient of 
Neameynenvwemaese surtace “Reéaghe “them “observed may be a 
combination of several factors of model dynamics. As 
mentioned earlier, the structure and magnitude of the actual 
headland winds is still in question and under investigation. 
The fact that the headland winds were modeled to oscillate in 
jee S “Oman vicarly cycle May Be Voo simplified cémoarcemto 
observations that show “upwelling favorable events” occur on 
a periodic basis, lasting perhaps several days, rather than 
remaining constant and diminishing on a monthly basis. The 
result of the enhanced upwelling may be such that the colder 
temperatures along the coast produce lower dynamic sea surface 
heights thereby increasing both the horizontal pressure and 
temperature gradients. This in turn may be the cause of the 
stronger than observed géostrophic velocities simulated along 
the upwelling front. 

To test the possibility that the headland winds are 
responsible for the stronger coastal currents (and SSH 
gradients) a control experiment should be conducted without 
the headland winds to see if the overall intensity of 
upwelling is reduced sufficiently to produce more realistic 


values and behavior of sea surface height and temperature 
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while maintaining the proper structure. To further enhance 
the simulation, the surface wind stress forcing with Hellerman 
and Rosenstein (1983), mean monthly climatological ind faetas 
couldybesreplacedsgwith data, £rompavhi gh mesolimiontatmospmensc 
model that shewSe ae merc realistwemsspemieiemwend temporal 
yarwabilatyeanethe wanewstress field. Sihastwould Velinenare 
the need for an idealized headland wind enhancement if the 
atmospheric model can properly produce the headland wind 
structure, inciuding 1tsesyneptic andysecasonalvariabiian, . 

The high levels of eddy kinetic energy west of about 
126°W to the western boundary are most likely caused by 
artificial interactions with the western boundary. As stated 
earlier, this may be the result of an improper (incompletely 
open) boundary condition not allowing the free and efficient 
exchange of energy across the western boundary. The cycle of 
the bazotropic andesaroclinuwe condasmyens ae thesnorehemm and 
southern boundaries are less problematic because eddy 
propagation is zonal (westward). Nevertheless, as seen in the 
reduced area/coastline relative domain, results can be 
interpreted devoid of the western boundary iméluence. 

One solution to the western boundary problem would be to 
make an ever bigger domain to increase the size of the buffer 
zone and thereby lessen the effects of the western boundary 


influence. This however would require significantly increased 
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pmece>cimG and data storage resources, not a very appealing 
Opie ON« Cneeeraeaaesolubion would) be to “nest” thé DIECAST 
model into an ocean basin model, like that of Semtner and 
Chervin aC eor a larger DIECAS? domain; where fully 
Ineeractime, boum@ary .condit ins "ear be usear" This wottta™be 
ahie moStwappecaling optien amd™ likely tomwyield «the most 
satisfactory results. 

In spite of the above shortcomings, the simulation showed 
a Pronounced tendency fer ciismero and downward propagdatrom of 
eddy kinetic energy with time, and the rate of westward 
propagation is consistent with observations. The downward 
spreading of EKE offshore etre samulation (Pigure Se 
offers a possible interpretation of the westward decrease of 
EKE at the surface found in the recent observational analyses 
by Kelly et al. (1998) and Strub and James (1998). Tires 
hypothesized that coastal eddies develop via baroclinic 
instability in the upper ocean in the vicinity of the 
continental slope, and as the eddies develop and propagate 
westward they transfer EKE to the deep ocean in a process 
Similar to the occlusion process in atmospheric cyclones. 
Further study is clearly needed to test such a hypothesis. 

The conclusion of this study is that the model simulation 
UmmreZ ing “advanced “second orders physics can *product“results 


consistent with observations and can provide insight to 
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interpret such observations. The inclusion of theseued ame. 
dynamics provides a significant improvement in the “realism” 
of the simulation over that carried out by Akahoshi (1995), 
but there is considerable room for additional improvements 


especially with regard to wind forcing and the treatment of 


open boundaries. 
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